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■ We present the analysis of the muon events with all muon multiplicities collected during 21804 

hours of operation of the first LVD tower. The measured angular distribution of muon intensity 
has been converted to the 'depth - vertical intensity' relation in the depth range from 3 to 12 km 
w.e.. The analysis of this relation allowed to derive the power index, 7, of the primary all-nucleon 
spectrum: 7 = 2.78 ± 0.05. The 'depth - vertical intensity' relation has been converted to standard 
rock and the comparison with the data of other experiments has been done. We present also the 
derived vertical muon spectrum at sea level. 

PACS numbers: 13.85.T, 96.40. T 



I. INTRODUCTION 



During the last 30 years the cosmic-ray muon energy spectrum has been studied in many experiments using different 
methods. These methods can be combined into 3_groups: i) the direct measurements of the muon energy spectrum 
at the sea level using magnetic spectrometers |l|-|3| (the spectrum up to about 10 TeV was measured at the zenith 
angles near horizon); ii) the measurements of the energy spectrum of cascadesproduced by muons at shallow depth 
H-0] ; iii) the measurement of the depth-intensity curve deep underground |^- [l^ , |33] | . 

Since the spectrum of primary cosmic-ray nucleons has the power-law form with the power index 7, the spectrum 
of 7T— and K— mesons produced by primaries should have also the power-law form with the power index j Wt K- If 



1 



the scaling hypothesis is valid in the fragmentation region at high energies, the value of 771-, if is approximately equal 
to 7. The muon spectrum at the sea level has a more complex form due to the competition between the interaction 
and decay of their parents. Moreover, the flux of muons produced by pions and kaons strongly depends on the zenith 
angle because the interaction path length of mesons in the atmosphere varies with the zenith angle 9. The muon 
spectrum follows the power-law dependence only at high energies E^» E%~, E^, where E% r and Eg are the critical 
energies of pions and kaons in the atmosphere. For 9 = (vertical) this condition is quite satisfied at E^ > 2 TeV. 
The power index of the muon spectrum, 7^, in this case is more by 1 than the value of Jtt,K- 

Despite the numerous experiments which studied the muon spectrum, there are the discrepancies in the published 
results (values of 7, "fn,K or 7^). Most of the experiments carried out in the last 20 years gave the values of 7 (or 
1k,k) in the range 2.60 - 2.80. But the dispersion of the results is greater than the statistical and systematic errors 
published by the authors. Even the values of 7, obtained using one method (for example, by the measurement of the 
depth-intensity curve) have a large dispersion. Such discrepancy can be due to either the difference between the data 
themselves, or between the calculated muon intensities used to fit the data, or both. 

The muon intensities deep underground have been calculated by many authors using different methods (Monte- 
Carlo simulation of the muon transport, numerical solution of the kinetic equations etc., see, for example, ||-|||). 
However, the discrepancy between the muon survival probabilities and, hence, between the muon intensities, calculated 
by different authors, is quite large and can result in the significant discrepancy between the final values of 7. 

The LVD (Large Volume Detector), located underground, can measure the atmospheric muon intensities from 3000 
hg/cm 2 to 12000 hg/cm 2 and above (which correspond to the muon energies at the sea level from 1.5 TeV to 40 
TeV) at the zenith angles from 0° to 90°. This allows us to study the muon spectrum and their characteristics at the 
energies 1.5-40 TeV (which correspond to the energies of primaries of about 10-400 TeV). 

In a previous paper p"2] | we have presented our measurement of the muon depth-intensity curve and the evaluation 
of the power index of the meson spectrum in the atmosphere using the depth-intensity relation for single muon events. 
The muon survival probabilities, used to obtain the value of j^.k in Jl2[ , have been presented in pl[ |. They have 
been calculated using the muon interaction cross-sections from |22|-|2~4|] . After the publication of these results, new 
calculation of the cross-section of muon bremsstrahlung and of the corrections to the knock-on electron production 
cross-section have been done pq| . In the present analysis we have taken into account the corrections proposed in 
pEf and we have estimated the uncertainties of 7 (7 is assumed to be equal to J-r^k) due to the uncertainties of the 
cross-sections used to simulate the muon transport through the rock. The analysis is based on an increased statistics 
comparing with the previous publications and refers to events with all muon multiplicities. 

In Section 2 the detector and the procedure of data processing together with the conversion of the muon intensity 
to vertical are briefly described. In Section 3 the results of the analysis of the 'depth - vertical muon intensity' 
distribution (I^(x)) are presented. In Section 4 the 'depth - vertical intensity' relation in standard rock is compared 
with the data of other underground experiments. In Section 5 we present the derived muon spectrum at sea level. 
Section 6 contains the conclusions. 



II. LVD AND DATA PROCESSING 

The LVD (Large Volume Detector) is located in the underground Gran Sasso Laboratory at a minimal depth of 
about 3000 hg/cm 2 . The LVD will consist of 5 towers. The 1st tower is running since June, 1992, and the 2nd one - 
since June, 1994. The data presented here were collected with the 1st LVD tower during 21804 hours of live time. 

The 1st LVD tower contains 38 identical modules J28|. Each module consists of 8 scintillation counters and 4 
layers of limited streamer tubes (tracking detector) attached to the bottom and to one vertical side of the supporting 
structure. A detailed description of the detector was given in p8[ . One LVD tower has the dimensions of 13 x 6.3 x 12 

o 

nr. 

The LVD measures the atmospheric muon intensities from 3000 hg/cm 2 to more than 12000 hg/cm 2 (which corre- 
spond to the median muon energies at the sea level from 1.5 TeV to 40 TeV) at the zenith angles from 0° to 90° (on 
the average, the larger depths correspond to higher zenith angles). 

We have used in the analysis the muon events with all multiplicities, as well as the sample of single muons. Our 
basic results have been obtained with all muon sample. This sample contains about 2 millions of reconstructed muon 
tracks. 

The acceptances for each angular bin have been calculated using the simulation of muons passing through LVD 
taking into account muon interactions with the detector materials and the detector response. The acceptances for 
both single and multiple muons were assumed to be the same. 

As a result of the data processing the angular distribution of the number of detected muons N^((f>, cos 9) has been 
obtained. The angular bin width 1° x 0.01 has been used. The analysis refers to the angular bins for which the 
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efficiency of the miion detection and track reconstruction is greater than 0.03. We have excluded from the analysis 
the angular bins with a large variation of depth. 

The measured A/^^, cos ^-distribution has been converted to the 'depth - vertical muon intensity' relation, I^{x), 
using the formula: 

, > _ £ t j N^Xmj^j, COS0j)) ■ cos9** 

" ~ (A-{x m {4>j > cos 6i))e{x m ((f>j , cos 6i)) • % • T) U 

where the summing up has been done over all angular bins (</>j, cos 9i) contributing to the depth x m ; A(x m (4>j, cos 6i)) 
is the cross-section of the detector in the plane perpendicular to the muon track at the angles (<pj, cos 9i); 
e(x m ((j)j ;, cos Qi)) is the efficiency of muon detection and reconstruction; f2y is the solid angle for the angular bin; 
T is the live time; and cos#** = I^(x m , cos = 1)/ I°(x m , cos 8i) is the ratio of predicted muon intensity at cos (9 = 1 
to that at cos 9i . To obtain the values of cos 6** the eq. (||) (see below) with the previously estimated parameters of 



the muon spectrum at sea level fl2|]21|j and the calculated survival probabilities have been used. Actually, the values 
of cos 9** do not depend much on the parameters of the muon spectrum at sea level (normalization factor A and power 
index 7) for any reasonable values of 7. The factor cos 9** is different from the simple cos 9 - law used to convert the 
muon intensities to vertical in Jl^ ]. 

In the calculations of cos#** the contribution of the prompt muons from charmed particle decay has been neglected. 



As it is shown in 26 27]] the ratio of prompt muons to pions according to LVD data does not exceed 2 • 10~ 3 at 95% 
confidence level. This means that at median depths (no more than 6-7 km w.e.) the fraction of prompt muons with 
respect to conventional muons does not exceed 10% at vertical. And there is also a little probability that this fraction 
is more than 10% at large depths (8-10 km w.e.). 

For depth - intensity relation the bin width of 200 m w.e. has been chosen. For depths more than 9 km w.e. the 
bin width was increased to 500 m w.e. to increase the statistics for each bin. The conversion of muon intensity to the 
middle points of each depth bin has been done using the formula: 



-^ Xi ) = i^( Xm )3j^L ( 2 ) 

J-^yXm) 

where I™(x m ) and I^(x m ) are the measured and calculated muon intensities at the weighted average depth x m 
which corresponds to the depth bin with the middle value of Xi\ I™(xi) and I^(xi) are the derived and calculated 
muon intensities at the depth Xi which is the middle point of the depth bin. The values of x m have been obtained 
by averaging the depths for all angular bins contributing to the given depth bin with a weight equal to the detected 
number of muons. To calculate the muon intensities at x m and Xi we have used again eq. (Q) with the previously 
estimated parameters of muon spectrum Jl2|,^lJ and the simulated muon survival probabilities. 

Since the width of depth bins is quite small (200 m w.e. for depth bins with high statistics) and the number of 
angular bins contributing to each depth bin is quite large (several hundreds), the conversion factor does not exceed 
10%. 



III. 'DEPTH - VERTICAL INTENSITY' RELATION IN GRAN SASSO ROCK 

'Depth - vertical muon intensity' relation derived as it is described in the previous section has been fitted with the 
calculated distributions with 2 free parameters of the muon spectrum at sea level: normalization constant, A, and the 
power index of the primary all-nucleon spectrum, 7. To calculate the muon intensities underground we have used the 
formula: 

r( n\ I™ v(T? ,dI^{E^cos9l 

I fl (x,cos9) = P{E^,x) dEf^o (3) 

JO Ctii/i0 

where P(E^q,x) is the survival probability of muon with an initial energy E^q at sea level to reach the depth x; 
and dI ^o(E^o,cos8) . g ^ e muon spectrum at sea level at zenith angle 9. This spectrum has been taken according to 



dl^Eyg, cos 9) = ^ Q U ^ 



dE^o 



7/0 



0.054 \ 
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where the values of cosO have been substituted by cosQ* which have been taken from pofl . In |30| cosO* = 
E^ K {cos9 = l)/E% r K (cos9), where E% r K are the critical energies of pions and kaons. This formula has been ob- 
tained under a simple assumption of scaling in the high-energy hadron- nucleus interactions. Under this assumption 
the power index of primary spectrum, 7, is equal to that of meson (pion + kaon) spectrum, Yt^k. To fit the 'depth - 
vertical intensity' relation measured by LVD we have put cos 9 = cos 9* = 1, however, to convert the muon intensity 
to vertical (as it was described in the previous section) we have used the values of cos9* from p0| . 

We have used the muon survival probabilities P(E^o,x) calculated with the muon cross-sections from p3|-p5| and 
with the account of stochasticity of all processes of muon interaction with matter. They differ from those presented 
in |2lj] and used also in |L2| (see also references therein) because of the new corrections for muon bremsstrahlung and 
knock-on electron production cross-sections proposed in |25) (see also the discussion in [^l] ) . 

The measured 'depth - intensity' curve has been fitted with the calculated function (see eq. (||)) with two free 
parameters: additional normalization constant, A, and the power index of primary all-nucleon spectrum, 7. As 
a result of the fitting procedure the following values of the free parameters have been obtained: A = 1.95 ± 0.31, 
7 = 2.78±0.02. These values are in good agreement with the results of the analysis of the depth - angular distributions 
p6jp7| . We note that the energy in eq. (Q) is expressed in GeV and the intensity is expressed in cm~ 2 s _1 sr _1 . 
The errors of the parameters include both statistical and systematic uncertainties. The latter one takes into account 
possible uncertainties in the depth, rock composition, density etc., but does not take into account the uncertainty 
in the cross-sections used to simulate the muon transport through the rock. If we add the uncertainty in the muon 
interaction cross-sections, the error of 7 will increase from 0.02 to 0.05 and the error of A to 1.0 (for the discussion 
about the uncertainty due to different cross-sections see pl|). Similar analysis performed for single muons reveals 
almost the same value of power index while the absolute intensity is 10% smaller: A — 1.65 ± 0.30, 7 = 2.77 ± 0.02. 
We note that the estimates of the parameters A and 7 are strongly correlated. The larger the value of 7 is, the larger 
the normalization factor A should be. The 'depth - vertical muon intesity' relation is shown in Figure 1 for all muon 
sample together with the best fit. The muon intensities are presented also in Table J (column 2). 

If the formula from Q is used for the muon spectrum at sea level instead of eq. (Q), the best fit values of 7 will be 
decreased by 0.04-0.05 and will be in agreement with the previously published values for single muons |l2|j2l|| analysed 
using the formula from J3^| . 

The value of 7 obtained with LVD data is in reasonable agreement with the results of many other surface and 
underground experiments (see, for example, |p|-P,p|,[7 pC| Jll| ] ) . 



IV. 'DEPTH - VERTICAL MUON INTENSITY' RELATION IN STANDARD ROCK 



The simulations carried out for Gran Sasso and standard rocks allow us to obtain the formula for the conversion of 
the depth in Gran Sasso rock, x gs , to that in standard rock, x s t ■ This was done by comparing the values of x s t and 
x gs for the same muon intensity: I^{x s t) — I^{x gs ). The muon intensities have been calculated with the value of 7 
which fit well the LVD data. The depth in standard rock can be evaluated from the depth in Gran Sasso rock using 
the formula: 

x st = -9.344 + 1.0063^ + 1.7835 ■ 10~ 6 :z; 2 s - 5.7146 • lO" 11 ^, (5) 

where the depth is measured in hg/cm 2 . This formula is valid for depth range 1-12 km w.e.. It has been used 
to convert the depth- muon intensity relation measured in Gran Sasso rock to that in standard rock to allow the 
comparison with the data of other experiments. 

The 'depth - vertical intensity' relation in the standard rock for all muon sample is presented in Figure 2. It can 
be fitted with a three parameter function: 

= A (^"exp^, (6) 
where A = (2.15 ± 0.08) • 10~ 6 cm" 2 s" 1 sr" 1 , x = (1155±^) hg/cm 2 , a = 1.93tg;^. 

The best fit function is also shown in Figure 2 by solid curve. The LVD data converted to the standard rock agree 
quite well with the best fit functions for the data of MACRO Jll| (dashed curve in Figure 2) and Frejus |Q (dotted 
curve which almost coincide with the solid curve) underground experiments, despite the difference in the formulae used 
for depth conversion. The LVD data do not contradict also to the function which fit the data of NUSEX experiment 
p0{ (dash-dotted curve). 

The muon intensities measured by LVD and converted to the standard rock are presented also in Table | (3rd 
column) . 
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V. MUON ENERGY SPECTRUM AT THE SEA LEVEL 



Using eq. (Q) with the estimates of the free parameters, the LVD data for Gran Sasso rock can be converted to the 
vertical muon spectrum at the sea level. The simple normalization procedure has been applied: 



l^x - 0, E m0 ) 



(7) 



where I™(x) and K(x) are the measured and calculated vertical muon intensities at the depth x in the Gran Sasso 
rock, I%q(x = 0,E m o) and I^ (x = 0,E m o) are the derived and calculated (using eq. (0)) differential vertical muon 
intensities at energy E m § at sea level. The median energies E m o(x) which determine the muon intensity at the depth 
x have been calculated using the equation: 

where dI ^ E ^° - is the muon energy spectrum at sea level. The values of I^ix — 0, E m0 ) derived from LVD data are 
presented in Table O. The derived (full circles) and calculated (middle solid curve) muon spectra at the sea level are 
shown in Figure 3 together with the data of MSU (diamonds) , ASD 0] (open circles) and the best fit of MACRO 
PI (dashed line). Upper and lower solid curves in Figure 3 represent the errors in the parameters and additional 
10% error in the absolute normalization of the muon flux. Muon intensities at sea level derived from LVD data are 
presented also in Table ||. 



VI. CONCLUSIONS 



The angular distribution of muon intensity measured by LVD has been converted to the 'depth - vertical muon 
intensity' relation. The analysis of this relation in the depth range 3000-12000 hg/cm 2 has been done and the 
parameters of the muon spectrum at sea level have been obtained: A = 1.9 ± 1.0, 7 = 2.78 ± 0.05. The errors include 
both statistical and systematic errors with the systematic error due to the uncertainty of the muon interaction cross- 
sections dominating. Similar analysis performed for single muon events revealed almost the same value of power index, 
while the absolute intensity is 10% smaller. The 'depth - intensity' relation has been converted to the standard rock 
and fitted with a three-parameter function. This relation agrees well with the data of other underground experiments. 
Using the measured 'depth - intensity' curve and the estimated parameters of the muon spectrum at sea level we 
have derived the vertical muon energy spectrum at sea level. 
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3 4 5 6 7 8 9 10 11 12 

depth (km w.e. of Gran Sasso rock) 

FIG. 1. 'Depth - vertical muon intensity' curve in Gran Sasso rock measured by LVD together with the best fit using eq. 
) with the parameters: 7 = 2.78 and A — 1.95. 
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3 4 5 6 7 8 9 10 11 12 

depth (km w.e. of standard rock) 



FIG. 2. 'Depth - vertical muon intensity' curve in standard rock. LVD data are presented together with the best fit using 
three-parameter function (see eq. (^)) and the best fits to the data of other experiments: MACRO [11] (dashed curve), Frejus 
[33] (dotted curve) and NUSEX [10] (dash-dotted curve). 
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FIG. 3. Vertical muon energy spectra at sea level: • - LVD data; o - ASD data [4]; o - MSU data [7]; solid curves - LVD 
best fit together with the curves representing the errors of the parameters and of the absolute flux normalization; dashed curve 
- MACRO best fit [11]. 
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TABLE I. Vertical muon intensities measured by LVD vs depth in the Gran Sasso and standard rocks. Errors include both 



statistical and systematic uncertainties. 



x, km.w.e. 


i M , (cm ssr) , G.b. rock 


i M , (cm ssr) , st. rock 


3.1 


(2.09 ± 0.02) • 10~ 8 


(2.17 ±0.02) ■ 10~ 8 


3.3 


(1.58 ± 0.02) ■ 10~ 8 


(1.65 ± 0.02) ■ 10~ 8 


3.5 


(1.18 ± 0.01) ■ 10~ 8 


(1.23 ± 0.01) ■ 10~ 8 


3.7 


(8.82 ± 0.09) ■ 10~ 9 


(9.24 ± 0.09) ■ 10~ 9 


3.9 


(6.66 ± 0.07) ■ 10~ 9 


(7.04 ± 0.07) • 10~ 9 


4.1 


(5.08 ± 0.06) ■ 10~ 9 


(5.37 ±0.06) ■ 10~ 9 


4.3 


(3.87 ± 0.04) ■ 10~ 9 


(4.12 ± 0.04) ■ 10~ 9 


4.5 


(3.00 ± 0.03) ■ 10~ 9 


(3.19 ± 0.04) ■ 10~ 9 


4.7 


(2.34 ± 0.03) ■ 10~ 9 


(2.47 ± 0.03) ■ 10~ 9 


4.9 


(1.83 ± 0.02) ■ 10~ 9 


(1.99 ± 0.03) • 10~ 9 


5.1 


(1.37 ± 0.02) • 10~ 9 


(1.50 ± 0.02) ■ 10~ 9 


5.3 


(1.07 ± 0.02) ■ 10~ 9 


(1.16 ± 0.02) • 10~ 9 


5.5 


(8.07 ± 0.13) ■ 10~ 10 


(8.99 ± 0.14) • 10~ 10 


5.7 


(6.42 ± 0.11) • 10~ 10 


(6.79 ± 0.12) • 10 -10 


5.9 


(5.07 ± 0.09) • 10~ 10 


(5.52 ± 0.10) ■ 10~ 10 


6.1 


(3.94 ± 0.08) • 10~ 10 


(4.34 ± 0.09) ■ 10 -10 


6.3 


(3.11 ± 0.07) • 10~ 10 


(3.53 ± 0.08) • 10~ 10 


6.5 


(2.48 ± 0.06) ■ 10~ 10 


(2.68 ± 0.07) • 10~ 10 


6.7 


(2.03 ± 0.07) • 10~ 10 


(2.27 ± 0.07) • 10~ 10 


6.9 


(1.66 ± 0.06) • 10~ 10 


(1.88 ± 0.07) • 10~ 10 


7.1 


(1.25 ± 0.05) ■ 10~ 10 


(1.40 ± 0.05) • 10 -10 


7.3 


(l.OO ± 0.04) ■ 10 


(1.14 ± 0.05) ■ 10 


7.5 


(7.70 ± 0.34) • I0~ n 


(9.32 ± 0.39) ■ 10 -11 


7.7 


(6.27 ± 0.30) • 10~ n 


(6.85 ± 0.32) • 10~ n 


7.9 


(4.68 ± 0.25) • 10~ n 


(5.65 ± 0.29) • 10~ n 


8.1 


(4.27 ± 0.25) • 10~ n 


(4.45 ± 0.24) • 10~ n 


8.3 


(3.19 ± 0.25) • 10~ n 


(3.80 ± 0.26) • 10~ n 


8.5 


(2.67 ± 0.26) • 10~ n 


(2.97 ± 0.25) • 10 -11 


8.7 


(2.26 ± 0.26) • 10~ n 


(2.51 ± 0.27) ■ 10~ n 


8.9 


(1.40 ±0.22) • 10~ n 


(2.11 ±0.26) • 10~ n 


9.25 


(1.23 ±0.17) ■ 10~ n 


(1.40 ±0.16) • 10~ n 


9.75 


(6.8 ± 1.3) • 10~ 12 


(7.9±1.2)-10~ 12 


10.25 


(3.0 ±0.8) • 10~ 12 


(4.5±1.0)-10^ 12 


10.75 


(2.5 ±0.7) • 10~ 12 


(2.7 ±0.8) ■ 10~ 12 


11.25 


(1.8 ±0.7) • 10~ 12 


(2.2 ±0.7) ■ 10~ 12 


11.75 


(1.2 ±0.7) • 10~ 12 


(1.4 ±0.8) ■ 10~ 12 
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TABLE II. Vertical muon spectrum at sea level derived from LVD data. Errors include both statistical and systematic 
uncertainties. 



E^o, GeV 


J M , cm s sr GeV 


1932 


(1.30 ± 0.01) • 10 -11 


2793 


(3.42 ± 0.03) ■ 10~ 12 


3909 


(9.72 ± 0.11) • 10~ 13 


5371 


(2.88 ± 0.04) • 10~ 13 


7246 


(<) 73 + 16V 10~ 14 


9685 


(3.24 ±0.07) ■ 10~ 14 


12770 


(1.15 ±0.04) ■ 10~ 14 


16750 


(4.74 ±0.32) • 10~ 15 


21980 


(1.62 ±0.22) • 10~ 15 


28580 


(6.7 ± 1.5) ■ 10~ 16 


42660 


(1.6 ±0.5) ■ 10~ 16 
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